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We present the Bayesian analysis of four different types of backreation models, which are based 
on the Buchert equations. In this approach, one considers a solution to the Einstein equations for 
a general matter distribution and then an average of various observable quantities is taken. Such 
an approach became of considerable interest when it was shown that it could lead to agreement 
with observations without resorting to dark energy. In this paper we compare the ACDM model 
and the backreation models with SNIa, BAO, and CMB data, and find that the former is favoured. 
However, the tested models were based on some particular assumptions about the relation between 
the average spatial curvature and the backreaction, as well as the relation between the curvature and 
curvature index. In this paper we modified the latter assumption, leaving the former unchanged. 
We find that, by varying the relation between the curvature and curvature index, we can obtain 
a better fit. Therefore, some further work is still needed - in particular the relation between the 
backreaction and the curvature should be revisited in order to fully determine the feasibility of the 
backreaction models to mimic dark energy. 

PACS numbers: 98.80.-k, 95.36.+X 



I. INTRODUCTION 

The Universe, as observed, is very inhomogenous on al- 
most all scales. However, in a standard approach to cos- 
mology, it is assumed that the Universe can be described 
by the homogeneous and isotropic Friedmann-Lemaitre- 
Robertson-Walker (FLRW) models. The FLRW models 
provide a remarkably precise description of cosmological 
observations, but to achieve this we need to pay a price - 
in order to obtain concordance with observations it must 
be assumed that the Universe is filled with an unknown 
substance called dark energy. However, this substance 
has never been observed directly and, since it has very 
unusual properties, some have begun to ask whether dark 
energy is real or if it is the description of the Universe 
which, requires the existence of such an exotic entity, that 
is invalid. 

While it is possible that our Universe is filled with dark 
energy, many alternatives have been proposed: brane- 
world cosmologies (see [l| for a review), f(R) cosmology 
(see for a review), application of inhomogeneous cos- 
mological models (for a review see [3]) and others. One 
of the recently proposed approaches is based on an aver- 
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aging framework. Such an approach is motivated by the 
fact that the Einstein equations are non-linear, which 
means that the solution of the Einstein equations for a 
homogeneous matter distribution is different from the av- 
eraged solution to the Einstein equations for a general 
matter distribution. In other words, the evolution of the 
homogeneous model might be slightly different from the 
evolution of an inhomogeneous Universe, even though in- 
homogeneities in the Universe, when averaged over a suf- 
ficiently large scale, might tend to be zero. The difference 
between the evolutions of homogeneous and inhomoge- 
neous models of the Universe is known as the backreac- 
tion effect. In this approach, one considers a solution to 
the Einstein equations for a general matter distribution 
and then an average of various observable quantities is 
taken. Under certain assumptions such an attempt leads 
to the Buchert equations The Buchert equations are 
very similar to the Friedmann equations, except for the 
backreaction term, which is in general nonvanishing if 
inhomogeneities are present. For a review on the back- 
reaction effect and the Buchert averaging schemem the 
reader is referred to 0-01 . Based on this^cheme, Larena 
et al. have recently proposed a model [8] in which the 
metric of the Universe at a given instant looks like the 
FLRW metric, but the evolution of the scale factor is 
governed by the Buchert equations. In this paper we 
aim to perform the Bayesian analysis of the cosmological 
observations within the models proposed in 
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II. BACKREACTION MODELS 

If the averaging procedure is applied to the Einstein 
equations then for irrotational, pressureless matter and 
3+1 ADM space-time fohation with a constant lapse and 
a vanishing shift vector, the following equations are ob- 
tained |J| 



3- = -47rG(p) 
a 



3^ = 87^G(p)-i(7^)-iQ, 



(1) 
(2) 
(3) 



where a dot (') denotes dt , {TV) is an average of the spacial 
Ricci scalar '^^^TZ, Q is the scalar of expansion, a is the 
shear scalar, p is the matter density, and ( ) is the volume 
average over the hypersurface of constant time: (A) = 
(/ cfix\/—h)^^ J d^x\/~hA. The scale factor a is defined 
cube root of the volume: 



1/3 



(4) 



where Vq is an initial volume. 

Equation ([1]) is compatible with © if the following 
integrability condition holds 



^dt (Qa6) -1- la, {{R)a') = 0. (5) 



Similarly, as in the FLRW models, the following param- 
eters can be introduced: 



d 8^G , „ _ W o - <2> 



(6) 



was first introduced by Paranjape and Singh though 
in their model k was constant. The motivation for k(t) 
comes from the fact that the averaged spatial curvature, 
if calculated at one instant, does not have to be the same 
as the averaged spatial curvature calculated at another 
instant. This is closely related to the fitting problem 
closely studied by Ellis and Stoeger [Tol|. In considering 
the fitting problem, it becomes apparent that a homo- 
geneous model fitted to inhomogeneous data can evolve 
quite differently from the real Universe. Therefore, if in- 
homogeneous model is averaged at one instant its FLRW 
parameters may be different for the FLRW parameters 
obtained after averaging at another instant. 

The Buchert equations do not form a closed system. 
To close these equations and thus to calculate the evolu- 
tion of the scale factor one has to introduce some further 
assumptions [Ij. One such assumption can be: {TZ) ^ Q 
Q. As seen from the integrability condition ([5]), this 
leads to 



(7^> = (7^)ia" and Q 



n + 6 



(9) 



where n is an arbitrary parameter. Now, the final step is 
to derive a relation between the average spacial curvature 
(TZ) and the curvature index k. In analogy to the FLRW 
models, the following relation can be proposed ^]: 



'■{TZ) 



k{z) 



(n + 6)(l-r>^)(l + z)-("+^) 
|(n-f 6)(l-f)™)| 



(10) 



The Hamiltonian constraints can then be written as: 



In Sec. IIIIB 2 1 we will modify the above assumption and 
test models with different relations between k and {TZ). 
Summarizing, the model considered in this paper is de- 
scribed by the metric ([5]) , but the evolution of the scale 
factor is governed by the Buchert equations. Employ- 
ing the assumptions ^ and (jlOp the evolution equations 
reduce to the following relation: 



TZ 



1. 



(7) 



Observe that il-ji + i7g can act like J^a- Moreover, if the 
dispersion of the expansion is large then Q can be large 
and as seen from ([3]), one can get acceleration (a > 0) 
without the need for dark energy. 

The template metric of the Universe - the metric which 
describes the averaged universe can be written as 



d.s^ 



1 - fc(t)r2 



dr2 - a(tfr' {dd^ 



• Mif') 



(8) 

A similar approach, i.e. to consider the template met- 
ric with a scale factor which evolves accordingly to the 
Buchert equations instead of the Friedmann equations. 



H = Ho^n^il + z)3 + (1 - a«)(l + z)- 



(11) 



This model is parametrized by two parameters: flm and 
n. The distance, using (|S]), can then be calculated by 
solving 



dr 

dz 



1 — fcr^ 



r!„(i + z)3 + (i-r!„0(i + ^)- 



(12) 



Larena et al. [8] tested this model with a likelihood 
analysis using the supernova and CMB data. They found 
that this model is in good agreement with observations. 
In the next section we will perform the Bayesian analysis 
of this model using the type la supernovae (SNIa) data, 
baryon acoustic oscillations (BAO) and the observation 
of the cosmic microwave background (CMB) radiation. 
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III. BAYESIAN ANALYSIS 

The model presented in the preceding section will 
be confronted with cosmological observations in the 
Bayesian framework using the CosmoNest code [llj [l^ 
which was adapted to our case. In the Bayes theory all 
that we know about the vector of parameters (9) of a 
given model (M) is contained in the posterior probabil- 
ity density function (PDF), which is given by 



P{e\D,M) 



P{D\e,M)P{e\M) 
P{D\M) '■ 



(13) 





where D denotes the set of data used in the analysis; 
P{D\9,M) is the likelihood function for a given model, 
which will henceforth be referred to as L{9)\ P{9\M) is 
the prior PDF, which enables us to include our previous 
knowledge (i.e. without information coming from the 
data D) about the parameters under consideration; and 
the last quantity P{D\A4) is the normalization constant, 
called the evidence (or marginal likelihood) , which is the 
most important quantity in the Bayesian framework of 
model comparison. The posterior PDF could be simply 
summarized in terms of a best fit value, which could be 
the posterior PDF mode (the most probable value of 9) 
or the mean of the marginal posterior PDF of a given 
parameter (0j), obtained by integration (I13p over the re- 
maining parameters. 



FIG. 1: Constraints from SN data (SNLS - upper panel. Union 
- lower panel). 68% and 95% contours are presented. 



are interested in a wide range of models, we choose the 
parameter n to vary over a wide range n € [—4, 4]. 

The results are presented in Fig. [TJ As seen, the higher 
rim J the higher n must be taken in order to obtain a 
satisfactory fit. Also, the Union data set puts tighter 
constraints on the allowed range of the model parameters, 
although it prefers higher values of r2„i and n. 



Parameters estimation 



Supernova data 



Firstly, we consider observations of Type-la supernova, 
which are taken from the Supernova Legacy Su rvey |J_5J 
(SNLS) and the Union Supernova Compilation [ig. Af- 
ter analytical marginalization over the Hq parameter the 
likelihood function is of the following form 



LsN{^m,n) cx exp 



N 



(14) 

where N is the number of data points [N — 115 for the 
SNLS sample and N = 307 for the Union sample), 0"^ 
is the observational error jl^l, Xi = jif^'^"'^ — nl""^ , with 
— M {mi is the apparent magnitude and M 
the absolute magnitude of SNIa) , is the observed distance 
moduh [25], /if = 5 logig Dl + 25, where Dl = cr{l + 
z) is the luminosity distance, and r is given by p2p . 

We assume flat prior PDFs for the model parameters 
over their entire assumed ranges. The allowed range of 
f2m, is the same as the one used to constrain the (p function 
of eq. dm) [l3|, i.e. G [0.2,0.6]. Finally, since we 



2. CMB data 

The second set of cosmological observations comprises 
measurement of the CMB angular power spectrum. Here, 
we present the analysis of the CMB data which is re- 
stricted only to fitting the positions of the first two peaks 
and trough. Such analysis therefore ignores the shape of 
the CMB power spectrum. One, therefore, should keep 
in mind that full CMB constraints are tighter(26|. 

If it is assumed that the early Universe (before and 
up to the last scattering instant) is well described by 
the FLRW model then the CMB power spectrum can be 
parametrized by [13] 



(15) 



where 1x^12, and are the positions of the first, second 
and third peaks, and 



la 



(16) 



where is the co-moving distance to the last scatter- 
ing surface and is the size of the sound horizon at 
the last scattering instant. The function 0™ describes 
the phase shift of the m-th peak and is mainly sensitive 
to pre-recombination physics. It depends on the baryon 



density (flb/i^), where h = _ffo/(100 Mpc km s~-^), on the 
ratio of the radiation to matter density at last scatter- 
ing [pr{Zf.) / Pm{z*) = 0.042(rim^^)~"^(2*/10^)], where 
is the recombination redshift, on the spectral index (n^), 
and on the density of the dark energy before recombina- 
tion. 

We fit the positions of the first and second peaks, and 
the first trough of the CMB power spectrum [18]. We 
assume that rig — \ and neglect the density of dark 
energy before recombination. r» = c/ HQr{z^,), r{z) is 



given by 




31500nbh^TcMB/2.7Ky 



R 

W^n^h^iTcMB/^JK)-^. We take Tcmb = 2.728 and 
emp loy the fitting formulae for provided in Ref. 



empii 



We use the position of the peaks and trough locations 
in the CMB power spectrum coming from the WMAP3 
measurements [18,]. Those values were obtained in an 
almost model independent manner. The power spectra 
was fitted by the model composed of a Gaussian peak, 
parabolic trough and a second parabolic peak. Such 
model contains eight independent parameters, including 
peaks and though positions. Constraints were obtained 
using Markov Chain Monte Carlo method. The values 
and uncertainties of these parameters are the maximum 
and the width of the posterior probability distribution 
function. Additionaly the uncertainties include calibra- 
tion uncertainty and cosmic variance. The likelihood 
function could be therefore written in the following form 



LcMB{^m,n,flbh'^,h) cx 



exp 



'h - 220.8 
~0?7 



/3/2- 412.4 
1.9 



530.9' 



3.8 



(17) 



We assume flat prior PDFs for the model parameters 
in the ranges (case 1): D.m G [0.2,0.6], n € [-4,4], h e 



[0.64,0.80] (using information from HST [201]), ilbh^ e 
[0.0203,0.0223] (using information from BEN [21]). 

The 68% and 95% contours of the posterior PDFs (after 
marginalization over h and Q,bh?), obtained in the analy- 
sis with the CMB data as well as in the analysis with the 
joint data set (SNIa+CMB) (here the likelihood function 
has the following form L = LsnLcmb) are presented 
in Fig. [5] (left panels). As seen both data sets prefer 
n > — 1. However, unlike SNIa, CMB data prefers lower 
values of r2„i, but still the best fits are consistent with 
each other. 

The marginal posterior PDFs are presented in Figs. [3] 
and m (black lines). As seen, preferred values of 
are higher than inferred within the standard cosmologi- 
cal model. The means together with the 68% Bayesian 
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FIG. 2: Left panels: constraints from SN (blue), CMB (green) 
and SN+CMB (black) data (SNLS+CMB - upper panel, 
Union+CMB - lower panel). Right panels: comparison of 
constraints from SN+CMB data (SNLS+CMB - upper panel, 
Union+CMB - lower panel) with different prior assumptions 
on the parameters h and Q,bh^\ case 1 (black), case 2 (red), 
case 3 (green), case 4 (blue), case 5 (orange). 



confidence intervals (credible intervals) and the poste- 
rior modes are presented in Table U and Table for 
SNLS+CMB and Union+CMB respectively. An inter- 
esting observation is that the mean of n is around 1. 
When n = 0, the H{z) does not depend on the spatial 
curvature. Hence, if only expansion history is considered, 
in such a case the backreaction effect behaves as the cos- 
mological constant. When n — 2, the H{z) depends on 
the spatial curvature as in FLRW models. 

It is interesting to see how constraints on the VL„i and n 
are changed after changing the prior information on the 
h and Q,bh^ parameters. We consider four different situ- 
ations: case 2 - we fix the value of the Hubble parameter 
to the best fit value obtained in |20|], i.e. h = 0.72; case 3 
- we fix the value of Vtbh? to the best fit value obtained 
in [2l|, i.e. ^.^h? = 0.0213; case 4 - we fix the values 
of both parameters {h = 0.72 and ^bK^ = 0.0213). We 
observe that, CMB and SNIa data considered together 
prefer larger values of Q.bh? than expected from BBN 
constraints. Because of this we additionally check how 
the results are changed after allowing the Vlbh^ parame- 
ter to take larger values, i.e. expand the prior range to 
{^bh?)* € [0.020,0.025] (case 5). Constraints from SNIa 
and CMB data on the f2„i and n parameters for the cases 
described above are shown in the right panel of Figure [5] 
The largest changes occur when the Hubble parameter is 
fixed (blue and red contours). The values of VLm and n 
are shifted to larger values and the constraints become 
tighter (especially for the matter density). On the other 
hand, fixing the value of Q.bh? does not give substantial 



5 






FIG. 3: Marginal posterior PDFs for model parameters 
(SNLS+CMB) with different prior assumption on h and Q^h^ 
parameters: case f (black), case 2 (red), case 3 (green), case 
4 (blue), case 5 (orange). 



lines), but expanding the allowed range of flft/i^ changes 
it slightly (orange lines) . Fixing the value of h does not 
influence the posterior PDF of ilbh^ (red lines). SNIa 
and CMB data considered together prefer larger values 
of ilbh^, even when an extension of the allowed range 
is considered (orange lines). As one can see, changes in 
the constraints on the flm and n parameters are more 
prominent when the SNLS data set is considered. The 
means together with the 68% Bayesian confidence in- 
tervals (credible intervals), and the posterior modes are 
presented in Table U and Table M for SNLS+CMB and 
Union+CMB respectively. 



3. BAO data 

In addition to the geometric measurements described 
above, we study constraints obtained from the measured 
dilation scale of the BAO in the redshift space power- 
spectrum of 46,748 luminous red galaxies (LRG) from 
the Sloan Digital Sky Survey (SDSS). The dilation scale 
is defined as 



-I 1/3 



Hiz) 



(18) 






FIG. 4: Marginal posterior PDFs for model parameters 
(Union+CMB) with different prior assumption on h and Qth^ 
parameters: case 1 (black), case 2 (red), case 3 (green), case 
4 (blue), case 5 (orange). 



changes in the constraints (green and blue contours). We 
can see that expanding the allowed range shifts the best 
fit values of f2„j and n upwards but does not improve 
the constraints (orange contours). This can also be seen 
on the marginal posterior PDF plots. One can addition- 
ally conclude that fixing the value of Qbh^ does not have 
any influence on the marginal posterior PDF for h (green 



where Dj^ is the co-moving angular diameter distance 
and H{z) is the Hubble parameter as afunction of red- 
shift. The measured value of the dilation scale at z — 0.35 
is 1370 ± 64 Mpc. It should be noted that the value of 
1370 ± 64 Mpc was obtained within the framework of 
linear perturbations imposed on the homogeneous FLRW 
background. Instead, such an analysis should be carried 
out within the framework of the model considered in this 
paper. Otherwise, we should be aware of possible sys- 
tematical errors. When the geometry of the spacetime 
is not FLRW the possible sources of errors are: (1) the 
sound horizon can be distorted and can be of a different 
size in parallel and perpendicular directions; (2) the ex- 
pansion rate can be different with respect to parallel and 
perpendicular directions; (3) the redshift distortions, if 
analysed within the inhomogeneous model, might lead to 
estimates different from those received within the stan- 
dard approach; (4) another source of error comes from 
the fact that in their analysis Eisenstein et al. converted 
the redshifts of LRG galaxies to distances assuming the 
ACDM model. Despite these uncertainties we proceed 
with the analysis to see how the BAO data can possibly 
constrain the data. As seen from Fig. [SJ the measure- 
ments of the dilation scale at z = 0.35 do not put tight 
constraints on the parameters of the model. At higher 
redshifts the constraints will be tighter, and thus an anal- 
ysis of the BAO within the backreaction model will be 
required when future observational data is available. The 
likelihood function for the BAO has the following form 



LBAoi^m, n, h) cx exp 



thcor /-jobs \ 2 
V ~ ) 



(19) 
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FIG. 5: Left panels: constraints from SNIa (blue), CMB 
(green), BAG (red) and SNIa+CMB+BAG (black) data sets 
(SNLS+CMB+BAG - upper panel, Union+CMB+BAG - 
lower panel). Right panels: comparison of constrains from 
SN+CMB+BAG data (SNLS+CMB+BAG - upper panel, 
Union+CMB+BAG - lower panel) with different prior as- 
sumption on h and flt,h^ parameters: case 1 (black), case 
2 (red), case 3 (green), case 4 (blue), case 5 (orange). 



We assume flat prior PDFs for the parameters within 
the ranges described above for case 1. The 68% and 
95% contours of the posterior PDF (marginaUzed over h 
and flbh^ ) for the joint constraints from the supernovae, 
CMB and BAO data (with the likehhood function of the 
fohowing form L — LsnLcmbLbao) ^-re presented in 
Fig. [S] (black lines) . 

The marginal posterior PDFs are presented in Figs. [S] 
and [7] (black lines). In comparison with the preceding 
section, the PDFs are similar except for n, which is now 
shifted to higher values. The means of the marginal pos- 
terior PDFs together with the 68% credible intervals are 
presented in TableUand TablelTlfor SNLS+CMB+BAO 
and Union+CMB+BAO respectively. For example, n 



+0.80 
-0.76 



based on Union+CMB+BAO is now 1.60 
parison to 1.02^0 77 inferred only from Union+CMB. We 
will come back to the relation between the curvature and 
backreaction effects in Sec. IIII B 21 We also check how 
the results are sensitive to changes in the prior informa- 
tion for the h and ^li,h^ parameters. We consider all cases 
described above. Observe that changes are most promi- 
nent when the value of h is fixed. The contours become 
narrower (red and blue contours on the right panels of 
Figure O , which is due to tigher constraints on the ftm 
parameter. The values of f2„i and n do not change sig- 
nificantly. When we expand the allowed prior range for 
the parameter ftbh^, the values of ^l^n and n are slightly 
shifted upwards, and the constraints on n become weaker 
(see also the posterior PDFs presented in Figure [5] and 



FIG. 6: Marginal posterior PDFs for model parameters 
(SNLS+CMB+BAG) with different prior assumptions on the 
parameters h and Q,i,h?: case 1 (black), case 2 (red), case 3 
(green), case 4 (blue), case 5 (orange). 



Figure[7]). Constraints on h are also changed in this case: 
its value is shifted upwards. 

An interesting point is that the PDF of Q.hh^ increases 
with Vlbh^ and does not seem to decrease in the consid- 
ered range. This implies that the maximum of the PDF 
is out of the range and might suggest a tension between 
SNIa+CMB, SNIa+CMB+BAO and BBN constraints on 



B. Models comparison 

1. ACDM vs the backreaction model 

In this section we present the comparison between the 
model considered above (which will be referred to as 
model 1) and the ACDM model (which will be referred 
to as model 0). In the Bayesian framework, models are 
compared not only by how well they fit the data, but 
also by their complexity (see [2^ for a review) . The best 
model from the set of models under consideration is the 
one with the greatest value of the probability in the light 
of data defined as 



PiM,\D) 



P{D\AU)P{M,) 
P{D) 



(20) 



where Mi is a model under consideration and P{Mi) is 
the prior probability of a model. If we have no foundation 
for favouring one model in the set of models under consid- 
eration over another, we usually assume the same value 
of the prior quantity for all models, i.e. P{Mi) — 1/K, 
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TABLE I: Means of the marginal posterior PDFs for the 
model parameters, together with the 68% credible intervals 
for the SNIa+CMB and SNIa+CMB+BAO data sets. The 



corresponding values for the posterior modes are presented in 
brackets. 





SNLS+CMB 


SNLS+CMB+BAO 


n 
h 


0.37t^-;j5(0.32) 
0.87+n7? (0.00) 
0.75lg°33(0.79) 
0.0219t°:S (0.0223) 


0.421-^(0.40) 
1.59+^-^^1-48) 
0.731^3(0-73) 
0.02 19+0 0003(0-02 23) 


n 
h 


0.43t-.-(0.43) 
L55+n«Q(1.48) 
0.72 

Q.0219l°°S;g33 (0.0222) 


0.431-^(0.43) 
1.70+0 70(1-68) 
0.72 

0.0219l°-.°r3 (0-0223) 


n 
h 


0.37t;j;^^(0.31) 
0.76tO:«^(-0.05) 
0.75l°;°|(0.80) 
0.0213 


0.411-^(0.40) 
1.46irri(l-17) 
0.721^2(0-73) 
0.0213 


n 
h 


0.421-.- (0.41) 
1.36tO;^t(1.15) 
0.72 
0.0213 


0.42l--(0.41) 
1.53l«:«^(1.29) 
0.72 
0.0213 


n 
h 


0.40+^-^^(0.34) 
1.27l°-l^(0.27) 
0.75l°;S]^(0.79) 

o.o243ti;!5E!E;5 (0.0248) 


0.44l--(0.45) 
2.0li;i°(2.00) 
0.73lH2(0-73) 
0.0243lo-.°EJ°56 (0.0245) 


TABLE IL Means of the marginal posterior PDFs for the 
model parameters together with the 68% credible intervals 
for the SNIa+CMB and SNIa+CMB+BAO data sets. The 
corresponding values for the posterior modes are presented in 
brackets. 




Union+CMB 


Union+CMB+BAO 


n 
h 


0.40t;j;«(0.36) 
1.02lo-^?(0.35) 
0.74lo°33(0.76) 
0.0219+00003(0.0223) 


0.43i;j-'(;^(0.42) 
l-60ir7°(l-14) 
0.72lEJ:°^(0.72) 
0.02 1910 0003(0-0223) 


n 
h 


0.42+^-^^(0.43) 
L4ll°-.f3(1.51) 
0.72 

0.0219«-«r3 (0.0223) 


0.431-^(0.42) 
1.57l«:«J(1.34) 
0.72 

0.02191-.S (0-0223) 


n 
h 


0.391^:15^(0.38) 
0.97^075(0.48) 
0.731°-°^ (0.74) 
0.0213 


0.42l--(0.42) 
1.48155-^^(1.40) 
0.72lH2(0-72) 
0.0213 


n 
h 


0.421-.- (0.41) 
1.25l°-.f2(1.07) 
0.72 
0.0213 


0.421--(0.41) 
1.391° j3 (1.23) 
0.72 
0.0213 


n 
h 


0.411-;;' (0.39) 
1.30t°-^t(0.83) 
0.75lg-.g^3(0.77) 
0.0244l|5-°^j;g (0.0249) 


0-44+0 03(0-44) 
1.92l«;«t(l-94) 
0.73lH2(0-74) 
0.0243l»-.S (0-0249) 




0.2 0.3 0.4 0.5 0.6 -4 -2 2 

"m n 



0.65 0.7 0.75 0.8 0.02 0.022 , 0.024 

FIG. 7: Marginal posterior PDFs for model parameters 
(Union+CMB+BAO) with different prior assumptions on the 
parameters h and Qi,h^'- case 1 (black), case 2 (red), case 3 
(green), case 4 (blue), case 5 (orange). 

■where K is the number of models. P{D) is the normaliza- 
tion constant. P{D\Mi) is called the marginal likelihood 
(or the evidence) and has the following form 

PiD\M,)^ j L{e)p{e\M,)Ae = E,. (21) 

It is convenient to consider the ratio of posterior prob- 
abilities for the models which we want to compare. If 
prior probabilities for those models are equal, then the 
posterior ratio reduces to the ratio of the evidences. This 
ratio is called the Bayes factor (By = Ei/Ej). The val- 
ues of Bij are interpreted as follows: < In Bij < 1 as 
inconclusive, 1 < In Bij < 2.5 as weak, 2.5 < In Bij < 5 
as moderate and Ini?^ > 5 as strong evidence in favour 
of a model indexed by i with respect to a model indexed 
by j. The evidence [eq. (piT) ] was calculated using the 
CosmoNest code. We generated five chains for each case 
to obtain the uncertainty in the computed value of the 
evidence. We performed our calculations using the dif- 
ferent prior assumptions for the parameters h and Vlbh?, 
which were described in the previous section. The values 
of the logarithm of the Bayes factor calculated for the 
ACDM model (model 0) vs model 1 (i.e. the model pre- 
sented in Sec. IIII A[) - i?oi - are presented in Table IIIII 
and in Figure [51 

The comparison in the light of the SNIa data does not 
give conclusive results - this data set has not enough in- 
formation to favour one model over another. After the 
inclusion of information coming from the CMB, there is 
strong (SNLS+CMB) and almost strong (Union+CMB) 
evidence in favour of the ACDM model over the inhomo- 
geneous one. When we include information coming from 
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TABLE III: Values of the logarithm of the Bayes factor In Bij 
j for different priors and different data sets. 



calculated for the model indexed by i and the model indexed by 



Case 


Data set 




InSoi 


InBai 


lnBo2 


lnBo3 


lnB3i 


In Bo4 


In B41 


1 


SNLS 






0.89 ±0.10 


-0.22 ±0.10 


1.11 ±0.10 


1.15 ±0.07 


-0.25 ±0.07 


1.06 ± 0.12 


-0.17 ±0.12 




Union 






0.48 ± 0.08 


-0.13 ±0.09 


0.61 ±0.10 


0.74 ± 0.09 


-0.26 ±0.08 


-0.69 ± 0.09 


-0.21 ±0.08 




SNLS4 


-CMB 




5.52 ±0.14 


3.78 ±0.27 


1.74 ±0.25 


1.73 ±0.21 


3.79 ±0.24 


1.87 ±0.31 


3.65 ±0.33 




Union- 


fCMB 




4.94 ±0.15 


3.28 ±0.13 


1.66 ±0.15 


1.59 ±0.22 


3.34 ±0.21 


1.73 ±0.22 


3.21 ±0.21 




SNLS4 


-CMB4 


-BAO 


4.77 ±0.12 


3.13 ±0.19 


1.64 ±0.17 


1.83 ±0.21 


2.93 ± 0.23 


1.78 ±0.16 


2.99 ±0.18 




Union- 


fCMB- 


fBAO 


4.23 ±0.08 


2.58 ±0.18 


1.65 ±0.17 


1.66 ±0.11 


2.57 ±0.12 


1.80 ±0.11 


2.43 ±0.12 



SNLS+CMB 
Union+CMB 
SNLS+CMB4 
Union+CMB- 



BAO 

-BAG 



4.52 ±0.20 
4.40 ±0.13 
4.09 ±0.10 
4.08 ±0.15 



3.17 ±0.24 
2.78 ±0.10 
2.53 ± 0.09 
2.52 ± 0.12 



1.35 ±0.16 
1.62 ±0.12 
1.56 ±0.09 
1.55 ±0.15 



1.50 ±0.10 
1.83 ±0.16 
1.73 ±0.10 
1.69 ±0.14 



3.02 ±0.21 
2.57 ±0.15 
2.36 ±0.10 
2.38 ± 0.12 



1.36 ±0.15 
1.60 ±0.12 
1.55 ±0.09 
1.69 ±0.13 



3.16 ±0.23 
2.80 ±0.10 
2.54 ±0.09 
2.38 ±0.10 



SNLS+CMB 
Union+CMB 
SNLS+CMB4 
Union+CMB- 



BAO 
-BAG 



6.41 + 0.18 
5.63 ±0.19 
5.67 ±0.06 
5.15 ±0.17 



4.62 ±0.17 
4.13 ± 0.23 
3.82 ± 0.20 
3.47 ±0.13 



1.78 ±0.20 
1.50 ±0.25 
1.84 ±0.20 
1.69 ±0.17 



1.77 ±0.27 
1.63 ±0.21 
1.89 ±0.11 
1.97 ±0.18 



4.64 ± 0.24 
4.00 ±0.19 
3.77 ±0.12 
3.18 ± 0.14 



1.82 ±0.15 
1.67 ±0.22 
2.02 ±0.08 
1.90 ±0.21 



4.59 ±0.11 
3.96 ±0.20 
3.65 ±0.09 
3.25 ±0.18 



SNLS+CMB 
Union+CMB 
SNLS+CMB4 
Union+CMB- 



BAO 

-BAG 



5.21 + 0.12 
4.98 ± 0.06 
4.90 ±0.09 
4.63 ±0.14 



3.99 ±0.13 
3.48 ± 0.12 
3.41 ± 0.16 
3.08 ±0.15 



1.22 ±0.10 
1.50 ±0.12 
1.49 ±0.17 
1.55 ±0.11 



1.39 ±0.11 
1.49 ±0.13 
1.60 ±0.11 
1.76 ±0.12 



3.82 ±0.13 
3.49 ±0.13 
3.30 ± 0.09 
2.87 ±0.16 



1.20 ±0.17 
1.38 ±0.22 
1.62 ±0.10 
1.76 ±0.18 



4.01 ±0.19 
3.59 ±0.22 
3.28 ±0.08 
2.88 ±0.21 



SNLS+CMB 
Union+CMB 
SNLS+CMB4 
Union+CMB- 



BAO 

-BAG 



2.85 + 0.12 
2.28 + 0.21 
1.88 + 0.21 
1.42 + 0.14 



1.26 + 0.27 
0.77 + 0.20 
0.89 ± 0.22 
0.73 ± 0.22 



1.59 ±0.27 
1.51 ±0.17 
0.99 ±0.15 
0.69 ±0.20 



1.53 ±0.12 
1.38 ± 0.20 
0.92 ±0.13 
0.58 ±0.21 



1.32 ±0.12 
0.90 ±0.22 
0.96 ± 0.20 
0.83 ±0.23 



1.77 ±0.14 
1.47 ±0.22 
0.94 ±0.14 
0.66 ±0.14 



1.08 ±0.14 
0.81 ±0.24 
0.93 ±0.21 
0.75 ±0.17 



SNLS 

SNLS+CMB 
Union+CMB 
SNLS+CMB+BAO 
Union+OMB+BAO 



FIG. 8: Values of the logarithm of the Bayes factor InBoi 
calculated for the ACDM model and the model indexed by i 
for five different prior cases and different data sets (SNLS - 
magenta, Union - cyan, SNLS+CMB - black, Union+CMB - 
red, SNLS+CMB+BAG - green, Union+CMB+BAG - blue). 



moderate for the SNLS+CMB data set and weak other- 
wise. When the value of h is fixed (case 2), the value of 
InB becomes smaller for all data sets and the evidence 
in favour of the ACDM model is moderate. On the other 
hand, fixing the value of ^ll,h^ (case 3) shifts the value 
of In B upwards and the evidence in favour the ACDM 
model becomes strong. Fixing both h and fl^h'^ (case 4) 
does not change the final conclusions. 

It is also interesting to see how strong is the influence 
of the prior information regarding the n parameter on 
the final results. We repeat our calculations for the case 
in which the prior range for the n parameter is changed, 
i.e. for n G [—2,3]. The values of the logarithm of the 
Bayes factor calculated with respect to the ACDM model 
for different data sets and different prior cases regarding 
the h and fli,h^ parameters are presented in Table HVl and 
in Figure [HI 

Observe that, when the prior range for n is restricted, 
the values of the logarithm of the Bayes factor be- 
come smaller, but the differences in the case in which 
n e [—4, 4] are small. Considering the comparison with 
the ACDM model, one finds that in general the final con- 
clusions do not change. 



the BAO, the values of the Bayes factor become smaller in 
both cases and the evidence in favour the ACDM model 
is moderate. 

Let us consider the influence of the prior information 
for the h and flhh^ parameters on the evidence. The most 
prominent change is when the range of f2(,/i^ is extended 
(case 5). The evidence in favour of the ACDM model is 



2. Relation between the average curvature and the 
curvature index 

In the preceding subsection we could see that the 
Bayesian method of model comparison prefers the ACDM 
model over the backreaction model. We should be aware 
that the backreaction model - model 1 - is based on 



9 



TABLE IV: Values of the logarithm of the Bayes factor In Boi calculated for the ACDM model and model 1 with assumption 
that n G [—2, 3] for different priors and different data sets. 



Data set 




case 1 


case 2 


case 3 


case 4 


case 5 


SNLS 




0.62 ±0.11 










Union 




0.17 ±0.14 










SNLS+CMB 




5.14 ±0.16 


4.12 ±0.13 


5.84 ±0.18 


4.79 ±0.10 


2.36 ±0.13 


Union+CMB 




4.48 ±0.17 


4.00 ±0.14 


5.13 ±0.17 


4.53 ±0.15 


1.86 ±0.18 


SNLS+CMB4 


-BAO 


4.32 ±0.15 


3.81 ±0.14 


5.32 ±0.05 


4.49 ± 0.18 


1.72 ± 0.14 


Union+CMB- 


fBAO 


3.80 ±0.11 


3.56 ±0.19 


4.61 ±0.19 


4.32 ± 0.12 


1.19 ± 0.18 



SNLS 
; SNLS+CMB 

Union+CMB 

SNLS+CMB+BAO 
■ Union+CMB+BAO 



sets, there is moderate evidence to favour the back- 
reaction model with k ^ over the model with 
relation PT!]) . Note that there is only weak evi- 
dence to favour the ACDM model over model 2. 
A change in the prior information regarding h and 
Slftft,^ does not change the final coirclusions for all 
considered cases, besides case 5. When the range 
for the parameter ^li,h^ is extended, the evidence 
against model 1 becomes smaller. 

model 3 



FIG. 9: Values of logarithm of the Bayes factor InSoi cal- 
culated for the ACDM model and Model 1 (with two dif- 
ferent prior assumptions on n parameter) for five different 
prior cases (regarding h and Qbh^ parameters) and differ- 
ent data sets (SNLS - magenta, Union - cyan, SNLS+CMB 
- black, Union+CMB - red, SNLS+CMB+BAO - green, 
Union+CMB+BAO - blue) 



the assumptions ^ and pO|) . If these assumptions are 
changed, it is possible to obtain a model with a better 
fit. 

Below we present models in which the assumption (jlOp 
is replaced with: 

• model 2 



k{z) = 0. 



(22) 



We emphasize that k — comes from a modifica- 
tion of the assumption pop only. It is not the same 
as assuming that (TZ) =0. As seen from (O the 
assumption of (7?.) — leads to Q ^ a^^, which 
means that the backreaction is strong in the early 
Universe and its value decreases with time. 

The results of the model comparison are presented 
in Table Hill and in Figure [S] As can be seen, both 
for the SNIa+CMB and SNIa+CMB+BAO data 



k{z) = -i(l + z)-("+2). 
P 



(23) 



We assume a fiat prior PDF for the additional pa- 
rameter in the range p g [0,100]. Observe that, 
model 2 fits the observations better than model 1. 
If indeed k = is favoured, then we should obtain 
that the best model with k given by (|23|) is tire one 
with p = 100. However, this is not the case and as 
seen from Fig. [TUl the marginal posterior PDF is 
almost flat over a very wide range - there is only 
a little difference between p = 100 and the best-fit 
value. This conclusion is coirfirmed bu observing 
the value of the logarithm of the Bayes factor: there 
is not enough informatioir to favour model 2 over 
model 3. The means and modes of the marginal 
posterior PDFs for tlie model parameters are pre- 
seirted in Table |Vl The values of the logarithm of 
the Bayes factor, calculated for the ACDM model 
and model 3 (In Bos), as well as for model 3 and 
model 1 (lni?3i) are compared in Table Hill (see also 
Figure [S]). One may conclude that there is weak ev- 
idence to favour the ACDM model over model 3 and 
moderate evidence in favour of model 3 over model 
1 (for the SNIa+CMB and SNIa+CMB+BAO data 
sets). The conclusion is chairged when the range for 
flbh"^ is extended: for the SNIa+CMB+BAO data 
set, the evidence in favour of the ACMD model be- 
comes inconclusive. 

model 4 



fc(z) =-(l + z)-("+2+™). 



(24) 
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SNLS+CMB SNLS+CMB+BAO 




20 40 60 80 100 20 40 60 80 100 

P P 



FIG. 10: Marginal posterior PDF for the parameter p of 
model 3 with different prior assumptions on the h and S7f,/i^ 
parameters: case 1 (black), case 2 (red), case 3 (green), case 
4 (blue), case 5 (orange). 



SNLS+CMB SNLS+CMB+BAO 




012345 012345 



Union+CWB Unior+CMB+BAO 




012345 012345 



FIG. 11: Marginal posterior PDF for the parameter m of 
model 4 with different prior assumptions on the h and Vlbh? 
parameters: case 1 (black), case 2 (red), case 3 (green), case 
4 (blue), case 5 (orange). 



We assume a flat prior PDF for the additional pa- 
rameter in the range m S [0, 5]. The marginal pos- 
terior PDF for the parameter m is presented in Fig. 
[TTl The means and modes of the marginal poste- 
rior PDF for the model parameters are presented 
in Table IVII The values of logarithm of the Bayes 
factor calculated for the ACDM model and model 4 
(lni?o4) as well as for model 4 and model 1 (lni?4i) 
are compared in Table Hill (see also Figure [S]). 

Observe that there is weak evidence to favour the 
ACDM model over model 4 for the all considered 
cases, besides case 5 and for the SNIa+CMB+BAO 
data sets, where the comparison is inconclusive. 
The differences in the values of the evidence be- 
tween model 4 and model 2 or model 3 are small, 
which leads to the conclusion that there is insuffi- 
cient information to favour the former. Model 4 is 
better than model 1 in light of the data used in this 
analysis, with moderate evidence against the latter. 
This conclusion is changed in case 5, in which the 
value of the Bayes factor calculated for those mod- 
els becomes smaller. 

Observe that except for the SNLS+CMB a change 
in the priors for h and Vlhh? does not change the 
shape of the posterior PDF for m, which is flat for 
m > 3. For the SNLS+CMB data, the posterior is 
wider, it is flat for m > 2 and becomes narrower 
when the value of h is fixed. When the range of 
is extended, the posterior PDF for m becomes 
wider in all cases. 

The above results are encouraging and motivate fur- 



ther study of the backreaction models. Especially, it is 
important to study assumptions other than ([9]) and ([TO)) . 
As seen, if only assumption ([TO)) was modified (models 2- 
4) then not only do we obtained a better fit, but also the 
value of Q.m realistically decreases compared to model 1. 



IV. CONCLUSIONS 

In this paper we presented a Bayesian analysis of the 
backreaction models. This work was motivated by the 
recently proposed model of the inhomogeneous alterna- 
tive to dark energy [8j]. In this approach the Universe is 
modeled by the Buchert equations, which describe the 
relations between the scale factor, the average spatial 
curvature, the average matter distribution, the average 
expansion and the shear. Larena et al. [8] showed that 
their model is consistent with supernova and CMB data. 
Here we included the BAO data and tested this model 
within the Bayesian approach (our model 1). 

Our analysis shows that the SNLS and CMB data 
alone strongly favour the ACDM model. With the Union 
sample and BAO data there is almost strong evidence 
(lni?2i > 4) to favour the ACDM model over our model 
1. However, if just the best- fit models are compared, then 
the (Union+CMB+BAO) for the ACDM and model 1 
are 320.96 and 325.36 respectively. If the distribution 
is assumed, then for 307 degrees of freedom in model 1 
the probability that this model is true in the light of data 
is 22.6%. In comparison, the ACDM model (308 degrees 
of freedom) gives 29.3%. Therefore, we can see that the 
best-fit model 1 fits observations almost as well as the 
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TABLE V: Mean of the marginal posterior PDF for the parameters of model 3, together with the 68% credible interval, for the 
SNIa+CMB and SNIa+CMB+BAO data sets. The corresponding values of posterior mode are presented in brackets. 





SNLS+CMB 


SNLS+CMB+BAO 


Union+CMB 


Union+CMB+BAO 


n 
h 

P 


0.24+-;;^ (0.20) 

-0.18l«i«(-0.37) 
0.76l°:»2(Q78) 

0.0216l«:«»^ (0.0222) 
53.18+3?-26 (75.80) 


0.28i;',;*(0.27) 

o-iolo::^f(-o.2o) 
0.741°:;;^ (0.74) 
o.o2i6iO:;]»;;^ (0.0221) 

54.25l^H? (63.71) 


0.2610:;;- (0.20) 

-0.13l°:3J(-0.53) 

„ no / „ \ 

0.75lO:ig(0.77) 
0.0216lO:»»»^ (0.0221) 

53.99135-19(50.23) 


o.3ii;;-oo7(o.3i) 

0.271°:™ (0.16) 

o.73i;;:;;^(o.73) 

0.02171»:000'^ (0.0223) 
53.87130 43(48.61) 


n 
h 

P 


0.30+^:^(0.28) 
-0.0ll°:«^(-0.35) 
0.72 

0.0216l[;'(;'(!'(!-^ (0.0220) 
56.14liEJi»(87.19) 


o.32i;;:i;^(o.32) 

0.30loj7"(-0-01) 
0.72 

0.0216lO:»™« (0.0222) 
54.35+igi^ (61.34) 


0.3210-0-5(0.34) 

0.18lg:°«(0.30) 
0.72 

0.0216+0-000^ (0.0223) 
55. 19liii;-ti (73.54) 


0.3310:0- (0.34) 
0.441^:^7(0.31) 
0.72 

o.o2i7i;i;;;;;i^ (0.0223) 
54.29+31-33(73.01) 


n 
h 

Qbh'^ 

P 


0.24+;',;^^3(0.20) 
-0.24l«:2S(-0.48) 
0.761°;°^ (0.77) 
0.0213 

53.09+31-l^q (89.47) 


o.27ro:t;"5(o.23) 

-0.0410*6 (-0.31) 
0.74lg:°3(0.76) 
0.0213 

55.331^°:^'^ (61.07) 


o.25i;;-;;:(o.2o) 

-0.23l0:^i(-0.57) 
0-7510:02(0-76) 
0.0213 

53.621^2:17(32.14) 


0.291°:- (0.24) 
0.0810:59 (-0.33) 
0.73l°:°!(0.75) 
0.0213 

55.03+3'i'-23 (84.36) 


n 

u 
n 

Clbh'^ 
P 


0.291^05(0.26) 
-0.19t°;g(-0.57) 

n 79 
u. / z 

0.0213 

56.95+1111(76.57) 


0.311^:^^5(0.26) 

o.oil°:r7(-o.54) 

n 79 
0.0213 

55.751^^-^^(78.65) 


0.301^-^^(0.27) 
-0.02l°-J«(-0.42) 
n 79 

0.0213 

55.84l|°:^o (91.07) 


0.3210.- (0.28) 
0.20lo-i|(-0.31) 
n 79 

0.0213 

55.37+;^;;-^^ (96.67) 


n 
h 

{nbhy 

P 


0.31+-;;^ (0.26) 
O-'iStlf 1(0.00) 
0.76lu;!3(0.78) 
0.0238i;^:^^J^ (0.0248) 
52.6ll^J;^^ (34.64) 


0.37+^-0^(0.36) 
1.28l}-:;«(1.14) 
0.74in^(0.75) 
0.02 40+0 0008(0-02 49) 
51-431^^:^^(76.04) 


0.34+0-0^(0.32) 
0.73lo-:|^(0.42) 
0.7410:03(0.76) 
0.0240+0 000? (0.0249) 
51-34lg:^^ (37.58) 


0.38l°:°o5(0.41) 

i.32i;;:9^(i.57) 

0-731°:;;^ (0.73) 
0.0241+^-;j;i;i^ (0.0248) 
51.44l32°2 (22.30) 


TABLE VI: Mean of the marginal posterior PDF for the parameters of model 4, together with the 68% credible interval, for 
the SNIa+CMB and SNIa+CMB+BAO data sets. The corresponding values of posterior mode are presented in brackets. 




SNLS+CMB 


SNLS+CMB+BAO 


Union+CMB 


Union+CMB+BAO 


n 
h 

Clbh^ 
m 


0.24+-"o3(0.20) 

-0.32l?,:29(-0.58) 
0.76+i;-"3(0.77) 

0.02161;;:™;;;^ (0.0220) 

3.00+}-^i(3.05) 


0.2810:0^(0.25) 

-0.06+;;-^ (-0.37) 

o.74i;;:»^(o.75) 
o.o2i6i»:;z;j (0.0221) 

3.20+}-i(3.18) 


0.26lOo:Og(0.23) 

-0.26lo-l°(-0.56) 
0.7410:0^(0.75) 
0.0216l°:°oo6 (0-0222) 
3.13lM^(3.79) 


0.3010:0^0.26) 

0-071°-;;° (-0.38) 
0-73i;i:;!!(0.74) 

o.o2i6i;!:;;;;;;t (0.0220) 

3.22+1-^^(3.61) 


n 
h 

Clbh^ 
m 


0.29+^:1^^(0.25) 
-0.24l°:g?(-0.64) 
0.72 

o.o2i5i;;:™« (0.0219) 

3.231} 1^(4.21) 


0.311^:00(0.29) 
0.04l°-^^(-0.33) 
0.72 

0.0216l":fHH,7 (0.0222) 
3.251^-1^(4.79) 


0.30lO:-(0.30) 
-0.10l0:g(-0.24) 
0.72 

0.021610:0007(0.0222) 
3.25+1-1!; (4.83) 


0.3210,- (0.31) 
0.17ir7^(-0.15) 
0.72 

0.0217i;i.;;°°t (0.0222) 
3.2211-1^(3.55) 


n 
h 

Q,bh? 
m 


o.24+i;:;j^3(o.2o) 

-0.38iE;-2,g(-0.55) 
0.76l°:|5^(0.76) 

0.0213 
3.00+1-33(4.69) 


0.2610:0^(0.22) 
-0.23+g-^|(-0.50) 
0.741°:°! (0.75) 

0.0213 
3.29+1-^^(4.50) 


o.25io:;;:(o.2i) 
-0.36+;;-;^^ (-0.64) 

0.7410:02(0.75) 

0.0213 
3.1711:^,(4.15) 


0.2810:0^(0.22) 
-0.1ll?,:;^t(-0.66) 

o.73i;;:;!!(o.74) 

0.0213 
3.30+1 19 (3.93) 


n 
h 

Clbh^ 
m 


0.28t;i:'^^(0.24) 
-0.39l°*^(-0.76) 
0.72 
0.0213 
3.26l}:f7(3.23) 


0.29l--(0.24) 
-0.16l°:°^(-0.73) 
0.72 
0.0213 
3.331}-}^ (3.66) 


0.291;;:;;- (0.25) 

-0.2710:^2 (-0.64) 
0.72 
0.0213 
3.3011:10(3.87) 


0.30l--(0.26) 
-0.05l°:«2(-0.52) 
0.72 
0.0213 
3.30+11^ (4.33) 


11 

h 

{fibh^r 

m 


o.3ii;;:;;?(o.25) 

0.37i!]:ro(-0.22) 
0.76li',-o3(0.77) 
0.0239l|;™'i'^ (0.0248) 
2.721^2(3-32) 


o.38i;;:o^(o.35) 

1.19ll:?i(0.70) 
0.741°:°^ (0.75) 
0.0240l°:°oo8 (0-0248) 
2.82l}-^«(4.10) 


0.34lo-o»(0.34) 
0.60lo-;|0(0.47) 
0.7410:02(0.75) 
0.02411°:°^!^^ (0.0249) 
2.7511:^^(3.55) 


0.3910:00(0.40) 

i.25i^,:;;j(i.43) 
0.731;;:;!" (0.74) 

0.0241+0-^;'^;^ (0.0249) 
2.7911:^1(3.23) 
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ACDM model. Still, there are other concerns regarding 
this model. For example, can the assumptions © and 
(fTO| be justified? In other words, how does the backreac- 
tion and the spatial curvature of the real Universe evolve, 
and does the relation Q ^ (TZ) for the real Universe hold? 
The most concerning is the value of Qm, which is quite 
large, 0.43to;o4. However, as it was shown in Sec. HUB 2| 
after the assumption (|10l) was modified, we were able to 
obtain a better fit and, in addition, the value of flm de- 
creased to 0.311q'q7 and O.SOlQ Qg (or even lower if the 
BAO data is excluded - see Tables fVl and I VII) for models 
4 and 5 respectively. This shows that a lot still needs to 
be done in the context of the backreaction models, es- 
pecially in the study of the relation between the average 
spatial curvature and the backreaction. 

We also investigated the sensitivity of the results to 
the prior assumptions on the h, n^h^ and n parameters. 
The most prominent changes are in the case in which the 
Hubble parameter is fixed and in the case in which the 
range of the parameter Jlb/i^ is extended. In the latter 
case, the evidence against model 1 with respect to the 
ACDM model is weak. 

Currently the ACDM model is preferred by the obser- 
vational data, but it is possible that, after the revision 



of assumptions © and (|10p we could obtain a more sat- 
isfactory results (see Table Hill and Fig. We should 
also remember that in these models of dark energy, the 
dark-energy-term appears as a consequence of inhomo- 
geneities that are present in the Universe. Therefore, 
within this class of models, the "decaying lambda term" 
takes on reveals a new and natural interpretation. 
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